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A mechanism of the quasi-one-dimensional (qld) superconductivity is investigated by 
applying the renormalization group techniques to the pairing interaction. With the ob- 
tained renormalized pairing interaction, the transition temperature and corresponding 
gap function are calculated by solving the linearized gap equation. For reasonable sets of 
parameters, of p-wave triplet pairing is higher than that of d-wave singlet pairing due 
to the one-dimensionality of interaction. These results can qualitatively explain the super- 
conducting properties of qld organic conductor (TMTSF)2PF6 and the ladder compound 
Sr2Cai2Cu2404i. 
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1. Introduction 

The nature of the attractive interaction and the type of Cooper pairmg in quasi-one- 
dimensional (qld) superconductivity have led to lively discussions in these decades. The firstly 
discovered and well-known qld superconductor is the Bechgaard salts (TMTSF)2PF6.^ Its 
superconducting (SC) phase is in proximity to the spin density wave (SDW) phase,2 so that it 
has been believed so far that the d-wave singlet pairing is realized due to the antiferromagnetic 
spin fluctuations.^"^ The recent experiments, however, are apparently against this expectation. 
The unsaturated behavior of the upper critical field Hc2 and the absence of the Knight shift 
suppression below the transition temperature on (TMTSF)2PF6 strongly suggests the triplet 
pairing. The Hc2 and the Knight shift data on the ladder material Sr2Cai2Cu2404i also 
suggest the triplet pairing. ^'^ 

In order to investigate the mechanism of the qld superconductivity, various theoretical 
attempts have been done. Previous studies on the qld Hubbard model have concluded the 
most favorable pairing symmetry is the d-wave singlet, which is reasonable in view of the 
unconventional superconductivity near the antiferromagnetic or SDW phase. On the other 
hand, some recent studies pointed out that triplet pairing can be realized assisted by charge 
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fluctuations. ^^"'^^ In a study based on third-order perturbation theory, which do not rely on a 
particular type of fluctuation but treat the moderate effect of electron-electron interaction, it 
was concluded that the d-wave singlet pairing is much more stable than the triplet one.^^ The 
quantum Monte Carlo study concluded that the /-wave triplet pairing can compete with (but 
cannot prevail over) the d-wave singlet one.^'' However, the validity of the existing theoretical 
methods are not so obvious in qld systems for the following reasons. 

(a) Perturbative corrections are not small at low energy and low temperature where logarith- 
mic singularities appear in each order of perturbation series. Therefore the higher-order 
corrections should also be taken into account. 



(b) Not only the specific fluctuations, such as spin or charge fluctuations, are developed, but 
also the vertex corrections have singularities. Therefore we should take into account each 
singular diagram on an equal footing. 

These are the effects of the one-dimensionality. 

The reason (a) suggests the low-order (third or even fourth order) perturbation theory 
is not enough, and (b) suggests that the random phase approximation (RPA) type approx- 
imation (including fluctuation exchange approximation), which emphasizes some particular 
fluctuations, such as spin fluctuations,^^ might lead to incorrect conclusion. 

The renormalization group (RG) techniques have been applied to one-dimensional (Id) 
systems as a powerful method which can overcome these difficulties.^^ In the conventional RG 
method for Id system (so-called g-ology) only two types of superconductivity, s-wave singlet 
and p-wave triplet, which are fully gapped on the Fermi surface, have been considered, and the 
superconductivity with gap having nodes on the Fermi surface have been out of considerations. 
It is, however, expected from the previous theoretical studies that there is a possibility that 
the gap has nodes on the Fermi surface in qld and even in purely Id systems. 

The purpose of this paper is to investigate the mechanism of the qld superconductivity by 
the newly developed RG method and to discuss the SC mechanism in repulsively interacting 
qld electron systems. This method can treat superconductivity with nodes, and can take into 
account the effects of the one-dimensionality. 

In §2 we present our theoretical method. The obtained gap functions and the transition 
temperatures for various situations are shown in §3. We give an interpretation of the results, 
and discuss the mechanism of the qld superconductivity in §4. Section 5 is devoted to con- 
clusions. We also present detailed momentum-dependence of the pairing interaction and a 
phenomenological discussion of a possibility of an odd-energy superconductivity in appendix 
A and B, respectively. 
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2. Theory 

2. 1 Outline 

The qld systems with a tiny interchain hopping (t^), i.e., nearly Id systems, are well- 
described by the Id RG. Suzumura and Fukuyama^*^ showed the t^-dependence of the tran- 
sition temperature, which exhibit the transition temperature of the SDW, TsdWj increases 
as TsDW oc for small and has a maximum at a certain t^, then turns down. (Here Q 
is a positive constant of the order of unity.) The superconductivity is realized above t^, but 
the Id RG cannot be applied for such relatively large tj^. So, we need a theoretical method 
which is appropriate for the qld system above t^, where the effects of the one-dimensionality 
is expected to be still remain. This is just what we shall present in this paper. 

The outline of the present method is as follows. First, we apply the RG technique to the 
irreducible part of the vertex, namely, the pairing interaction. With this renormalized pairing 
interaction, we solve the linearized gap equation 



in the weak-coupling formalism to determine the transition temperature and the corre- 
sponding gap function A(k) without any assumption on the functional form of the gap. 

Solving the linearized gap equation is equivalent to studying the Cooper instability with 
the sum of all matrix elements which are irreducible in the particle-particle channel in the 
sense that they cannot diagrammatically be divided into two parts by cutting two parallel 
electron lines. The pairing interaction y(k, k') thus has to be irreducible with respect to the 
particle-particle channel, otherwise the particle-particle channel is multiply counted when we 
solve the gap equation. We need to calculate not the full but the irreducible vertex, which 
cannot be obtained by the conventional RG method. 

2.2 Multiplicative RG technique 

The procedure of the multiplicative RG method is as follows. The multiplicative property 
of the Green function, vertex function and other physical quantities (e.g., response functions) 
allows one to improve perturbative results by solving the RG Lie equation. Assume that a 
quantity A{uj/EQ,gi) has a scaling property and satisfies 




(1) 






(3) 
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Fig. 1. Diagrams of the irreducible vertices P^n -^^ ^ up to -^j^ 2- The diagrams enclosed by [ ] have 
contributions i^-KVp)^^ \tl{uj jEo). The imaginary parts are neglected. 



where x = lo/Eq, T/Eq or vph/Eo, and the prime sign indicates the renormalized quantities. 
Notations for coupling constants here follows the conventional g-ology.^^ The scattering pro- 
cess with gi corresponds to backward scattering, whose momentum transfer is of the order 
of 2kF. The processes with g2 are forward scattering. The couplings with parallel (opposite) 
spins are indicated by the subscript || (_L). 

The first step of the RG method is to calculate the Green function and vertices perturba- 
tively. Using these perturbative forms in the right-hand side of the Lie equation and solving it, 
we obtain the renormalized solution for the invariant couplings. Once the renormalized invari- 
ant couplings are determined, we can also obtain other quantities (e.g., response functions) 
by solving the respective Lie equations. 

In the present RG method, the first step, where the renormalized invariant coupling is 
obtained, follows from the conventional g-ology.^^ The essential and specific point of our RG 
method is the second step, where we do not calculate the renormalized response functions, but 
the renormalized pairing interactions (which is two-particle irreducible in the particle-particle 
channel) . 

2.3 Irreducible vertices 

A relation between the irreducible part of the four-point vertex T* and the full vertex is 
given by the equation 

ri{uj) = T*{uj) + ij ^^T*{kF,-kF,p,-p;3uj/2,-u;/2,uj',uj -uj') 

'>iG{p,uj')G{—p,uj — u)')Ti{p, —p, —kp, kp'jOj' ,u! — w', a;/2, u;/2), 

(4) 
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where F's arc the fuh vertices and G is the fuh Green function. Suffices of F's have usual 
meanings. The irreducible parts are given, up to the lowest-order corrections, i.e., 0{gf), in 
the form 



^In 








Eq 


LO 


vri 


Eq 


~ 2" 



n±(^) = gi± + :^{gi\\9i±- 9i±92)\^n^- ^ j, (6) 

The corresponding diagrams are shown in Fig. 1. In order to consider not only these second- 
order contributions, but also the higher-order ones, we employ the RG technique. 

Introducing the dimensionless irreducible vertices F* as F* = giV*, the scaling differential 
equations (cf. eq. (3)) are obtained (up to one-loop order^^) as 

which are different from the scaling equation for the full vertices Fj.^^ These scaling equations 
(8)-(10) can be solved in a closed form in the case gi^^ = gij_ with use of the g-ology results, 
g'li^) = - {gi/'n:vF)lnx] and c/^(x) = 32 - gi/2 + g[{x) /2, as 

fl{x) = [l-(5i/7rt;i.)lnx]-'/V", (11) 

f*{x) = [l-{gi/7rvF)lnxf'x-''/^, (12) 

where a = {g2—gi/2)/TTVF. Near the fixed point, i.e., in the low-energy region, the logarithmic 
correction is less dominant than the power term, and is considered to be weakened when we 
apply higher order RG (e.g. two-loop RG). The renormalized irreducible vertex in the low- 
energy region is thus given by 

n ^ i^/Eo)-'', (13) 

F* {u/Eq)-''/\ (14) 

(The effect of the logarithmic correction will be discussed in §4.) Parenthetically, we note that 
the full vertices Fi^2 are renormalized as Ti{uj) 0,F2(i^) g2 — 5i/2 (constant). This is 
due to the fact that the zero-sound type and Cooper type singularity compensate each other. 
By contrast, the renormalized irreducible vertex F| 2 have the power-law singularity. This is 
understood as the lack of the compensation. (The diagrams which are effectively contained in 
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the renormalized T* are displayed in Fig. 7.) 

Using these renormaUzed irreducible vertices, we can construct the pairing interaction 
both for singlet and triplet channel. 

2.4 Pairing interactions for quasi- one- dimensional systems 
The pairing interactions are given by 

V^\ik-\-k') = -5,||(fc + fc') + 52||(A; + A;'), (15) 
Vlik-\-k') = -g,^{k + k')+g2±{k + k'), (16) 



for triplet, and 



V'{k-^k') = gi±{k + k')+g2±{k + k'), (17) 



for singlet. (See Appendix for details.) Using the results of the renormalization applied to 
the irreducible vertices, (13) and (14), the pairing interaction for 5111= gi± is given by 

With this pairing interaction, we can determine the SC transition temperature Tc and 
the gap function without any assumptions on the form of the gap functions by solving the 
linearized gap equation (1) numerically. Here, we take the qld dispersion 

^ (k) = —2t cos kx — 2t± cos ky — fi, (20) 

where t (t±) is the intrachain (interchain) transfer integral, and /j, is the chemical potential. 
Of course, when we solve the gap equation (1) in qld systems, we have to take a qld pairing 
interaction, i.e., V(k, k') must have dependence on kx, ky, k'^ and k'y. It is, however, naturally 
expected that the dominant character of the pairing interaction is determined by the Id pairing 
interaction. Following this idea, we take the quasi-one dimensionality into account as follows. 
For a qld dispersion ^{kx,ky) = vp{\kx\ — kp) — 2t± cos ky — 11, the logarithmic singularity in 
the Cooper and the zero-sound channels are cutoff as 

i^T^lnfeV (21) 
2TIVF \Eq ) ^ ' 

so that the qld Cooper and zero-sound channel may have a form 

2'itvf \ ko 

where q = \k -\- k' — 2kF\- Namely, the interchain hopping suppresses the singularity of the 
interaction (see Fig. 2). The renormalized irreducible vertices are also suppressed, and then 



^TT^^M ""^"^ (22) 
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-LO -0.5 0.0 0.5 LO 

Fig. 2. Model pairing interaction in qld system. The solid line is the qld interaction suppressed by 
the interchain hopping t}_. The dashed line indicates the Id interaction which has a power-law 
singularity at g = 0. 



we obtain the qld pairing interaction as follows: 



-a/2 



(23) 
(24) 

Although the pairing interaction (23) and (24) depends only on the momentum parallel to 
the chain, the quasi-one-dimensionality is taken into account. 

Using the pairing interaction, (23) and (24), and the qld dispersion (20) in (1), we calculate 
the transition temperature Tc and the gap function at Tc- We focus on the case of 1/4 filling 
as a model of (TMTSF)2X compounds, and take kp = 7r/4 (/cq = '^/^)- We take finer mesh 
near the Fermi surface, i.e., 1000 points in the region < < tt and the minimum size of 
the mesh is A/c™"^ = 6.0 x 10^^ (the energy resolution is 1.7 x 10^ ''t ~ 0.5mK). 

In the present framework, the gap equation gives finite Tc even for = 0. This is due to 
the mean field approximation for the SC instability, namely, the present method fully takes 
into account the zero-sound type diagrams, but partly does the Cooper type diagrams. For 
t± ~ 0, there is an interference between them generally, which leads Tc to zero.^*^ The region 
in which the interference is remarkable is T > Therefore the valid region of the present 
theory is T < ^ t. This region seems to be narrow, but the unknown and interesting 
phenomena happen there. The behavior for t± < T <^ t would smoothly continue to the 
g-ology with finite t±?^ 
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Fig. 3. Variety of the superconductive gap symmetry in qld systems, (left panels) The vertical line 
a.t k = kp = 7r/4 indicates the Fermi point in purely Id case {t± = 0). The right panels show the 
Fermi surface, where the solid (dashed) line indicates A{kx) > (< 0). 



3. Results 

The gap function of qld systems are roughly classified as shown in Fig. 3. It is true that, 
in qld systems, we cannot classify the gap symmetry by the spherical-wave basis such as 
s,p, d. We use such terminology to express how many times the gap changes sign on the Fermi 
surface, i.e., "s-wave" does not change its sign, "p-wave" changes twice, "d-wave" four times, 
and so on. Although the gap functions depend only on kx since the pairing interaction V^'^, 
does, the interchain hopping t_i slightly warps the Fermi surface, and the d- and /-wave gap 
thus have 4 line nodes on the Fermi surface. In other words, it seems sufficient for these four 
types of pairing symmetry to consider only fc^-dependence of the pairing interaction. 

The results of calculation are shown in Fig. 4 for various values of t_L and gi and for 
fixed g2 = O.Sttvf. The region we are interested in is the SDW (CDW) phase in the context 
of g-ology, i.e., < gi/g2 < 2.^^ In this parameter region, we obtain two types of pairing 
symmetries, p-wave triplet (pTS) and d-wave singlet (dSS). For t± = 0.02t, dSS is dominant 
for < gi/g2 < 1-0 and pTS is for 1.1 < 51/52 < 1-7. (Fig. 5) In the case of gi = g2, which 
corresponds to the Hubbard model, Tc of dSS is higher than that of pTS. When we increase 
t±, however, Tc of dSS is suppressed more drastically than that of pTS, so that pTS is the 
most stable for t± ^ 0.03t as shown in Fig. 6. Figure 4 and Fig. 6 tell us that dSS is "weaker" 
against tj_ than pTS, i.e., against the suppression of the singularity in the pairing interaction. 

In the present theoretical framework, the Id fluctuation is considered in the pairing in- 
teraction, but the SC transition is determined within the mean held approximation. When 
we consider the Id fluctuation effect upon the transition, it is expected that the T^s are 
suppressed as Tc ~ tj,^" as schematically shown by the dashed lines in Fig. 6. The SDW 
phase would overhang for t±/t ^ 0.1.^^ Considering the long-range nature of the Coulomb 
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Fig. 4. Interchain hopping t^/t - backward scattering 51/(72 - transition temperature Tc phase dia- 
gram for 52 = Q.^TTVp. The open circles indicate Tc of d-wave singlet superconductivity (dSS), the 
open triangles indicate that of p-wave triplet superconductivity (pTS). 
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Fig. 5. Superconducting transition temperatures as functions of gi/g2 for t±/t = 0.02 (top) and 
for t±/t = 0.1 (bottom). The circles indicate Tc of d-wave singlet superconductivity (dSS), the 
triangles indicate that of p- wave triplet superconductivity (pTS). 



interaction, we may take gi ^ g2, t±/t ~ 0.1 for the observed SC state in (TMTSF)2PF6. 
According to our calculation for t±/t = 0.1, the pTS arises for gi/g2 ^ 0.7. Even in such 
case, the field-induced pTS is realized by applying Zeeman magnetic field, while the ground 
state is dSS for 51/52 ~ 0.7, since the singlet pairing is suppressed due to the paramagnetic 
eff'ect while the triplet pairing is not. The signs of the triplet pairing in (TMTSF)2PF6 are 
mainly observed under the magnetic field. So, there is a possibility that the field-induced 
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Fig. 6. Interchain hopping t±_ dependence of Tc's in the case of gi = .92 = QSttvf, the Hubbard model 
case. The dashed hnes indicate the suppression due to the Id fluctuation. 



pTS is realized in (TMTSF)2PF6. The present results are also consistent with the experiments 
of Sr2Cai2Cu2404i, which shows coherence peak at and exponential decrease below in 
1/Ti suggesting the gap.^'^ (pTS is fully gapped for the present open Fermi surface.) 

4. Discussion 

Our result that the pTS instability is rather strong for parameters appropriate to the 
Hubbard model is different from other theoretical studies as described in §1. In order to 
clarify the difference between the present study and the previous ones, let us concentrate on 
the Hubbard model case {gi = 92 = U) for the moment. 

The diagrams of the pairing interaction for the Hubbard model on the basis of the pertur- 
bation expansion up to the quadratic order in logarithmic singularity are displayed in Fig. 7. 
Roughly speaking, negative promotes the p-wave triplet pairing, and positive promotes 
the /-wave. Likewise negative promotes the s-wave singlet pairing, and positive pro- 
motes the d-wave. (See Appendix .) On the basis of RPA-like approximation, which consider 
the diagrams enclosed by gray-solid lines in Fig. 7, we have 

yt'RPA^^) = U^Xoiq)>0, (25) 
V^'^'^'Hq) = 2U + U\oiq) + 2U'xliq)>0, (26) 

where q = k + k' and Xo('?) = — (l/2vrui7) In |(g — 2kF)/2kQ\ is the zero-sound or particle- 
hole channel, which corresponds to the bare spin susceptibility.^^ These form indicate that 
/-wave (d-wave) is stable for triplet (singlet) pairing, and the dSS is more favorable than /TS 
due to the large U^X^iq) term, which is consistent with the previous studies based on the 
antiferromagnetic spin fluctuation theory. ^^"^^ When we consider vertex corrections as well, 
we have 

V\q) = U\o{q)-2U'xliq) + U^XoiQ)<0, (27) 
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VHq) = 




Fig. 7. Diagrammatical expression of the effective pairing interaction for the Hubbard model on the 
basis of the second-order logarithmic perturbation. The bare spin susceptibility xo{<l) corresponds 
to the zero-sound channel, i.e., xo{q) — —{l/2'KVp)ln\{q — 2kp)/2ko\. {xoil) = Xo{q)/2ttvp) 
Diagrams enclosed with gray-solid line are included in the RPA-type approximation, and that 
with gray-dashed line are not included in the present one-loop approximation. 



V'{q) = 2U + U\o{q) + U'xo{Q)>0, (28) 

where Xo{q) = Xq{q) I'^'^vp- Since the bare spin susceptibility xo diverges at g = 2kF-, even 
if U is small, the relation < holds near q ~ 2k p- For the singlet pairing, V^{q) is 
still positive, but is smaller than V^'^^^{q), because the vertex corrections cancel out the 
large positive U^Xoil) term. Therefore, of the d-wave pairing is drastically suppressed 
by the vertex corrections, which is consistent with the third-order perturbation study. In 
this sense, the vertex corrections are important for the triplet pairing to prevail over the 
singlet pairing. The vertex corrections, moreover, can shift from /-wave to p-wave within 
the triplet pairing. In the case of the triplet pairing, V^{q) becomes negative, so that the 
pairing symmetry changes from /-wave to p-wave due to the most singular —2U'^x^{q) term. 
Even if we consider the suppression of Xo(q) due to t±_ or temperature, Xx){q) will remain 
enhanced compared to Xq{q) lid systems, i.e., the first and the third terms in eq. (27) 
cannot cancel out the second term, so that p-wave triplet pairing is considered to be robust. 
This consideration based on the second order logarithmic (third order in U) perturbation 
does not tell us that the vertex correction of 0(U'^) makes an essential contribution, but only 
tells us the singular vertex correction has a potential to change the pairing symmetry. This 
is the very Id effect. In order to clarify its influence, however, we have to examine higher- 
order term of the singular vertex corrections. The validity of the simple perturbation theory 
becomes more questionable due to the logarithmic singularity when the one-dimensionality is 
strong. The present results are obtained by considering both the vertex corrections and the 
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higher-order logarithmic singularities within the one- loop RG. The one- loop approximation is 
sufficient for the case > l^i^l, and thus the present method is appropriate to evaluate the 
higher-order leading terms. (Note that our calculation may overestimate this Xoil) term since 
one-loop approximation does not contain the next-leading diagrams enclosed by gray-dashed 
line in Fig. 7.) 

Incidentally, the logarithmic term in eq. (11) and (12) may not be negligible within the one- 
loop RG. In qld systems, the power-law singularity is suppressed as A(q) ~ {q'^ /'i+t'j_/vp)~"^'^ , 
which is comparable to the logarithmic singularity near q ~ 2t±/vF. If these logarithmic 
corrections are relevant, it suppresses the backward scattering F^ and enhances the forward 
scattering F2, namely, Tc of pTS is suppressed and that of dSS is enhanced. Though it is, 
this suppression (enhancement) for pTS (dSS) is considered to be small since the logarithmic 
corrections are generally irrelevant when we apply the higher-loop approximation. 

5. Conclusion 

We applied the multiplicative renormalization group method to the irreducible vertex as 
the pairing interaction and calculated the transition temperatures and the gap functions. 
The renormalized pairing interactions have the power-law singularity as Fl oc {q/2ko)~°', 
F2 oc (5/2/^0)^"^^, here a = (52 — gi/'^)/'^VF- Comparing their exponents, we may conclude 
that the irreducible vertex for backward scattering is much singular than the forward one. 
The solution of the gap equation shows that the p-wave triplet pairing can be realized in qld 
system {t±/t ~ 0.1, 51/52 ~ 1-0), prevailing over the d-wave singlet pairing. When we take into 
account the logarithmic correction to the irreducible vertices, Tc of the p-wave triplet state can 
be suppressed and that of the d-wave singlet state is enhanced. This logarithmic correction 
might be an artifact of the one-loop approximation, so we need further investigation whether 
this logarithmic correction is irrelevant or not. The two- loop approximation of the irreducible 
vertex, however, does not obey the scaling assumption, and thus we must take another route 
to settle this issue. Nonetheless, our new method allows us to recognize the one-dimensional 
effect which enhances the p-wave triplet instability in qld systems, which has not been noticed 
in the previous theoretical works. The obtained results may give interpretations of the spin- 
triplet behavior of qld materials (TMTSF)2PF6 and Sr2Cai2Cu2404i. 
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Fig. A-1. Two different diagrams of zero-sound type, (a) zero-sound channel contributing to the 
backward scattering; (b) that to the forward scattering. 



Appendix A: Pairing interactions for singlet and triplet pairing 

Here we mention the momentum dependence of the interaction. The first correction to the 
irreducible vertex is the zero-sound type diagram, which give contributions 

1 , vp\k ± k' — 2kp\ , , , 

-In— —. (A-1) 



27rvF Eq 

As seen in Fig. A-1, the zero-sound type contribution to the attractive interaction in the 
BCS Hamiltonian can be classified into two groups 

Xo{k - k')hl,^a)_j^,^,b_k„,aka, (A-2) 

XQ{k + k')a\,J)^_j^,^,b_k„,aka (A-3) 

Equation (A-2), displayed in Fig. A-1 (a), corresponds to the backward scattering, and eq. 
(A-3), displayed in Fig. A-1 (b), to the forward scattering. Therefore, we obtain the pairing 
interaction of the present Id model as follows: 

k,k' 

+ 9i±{k - k') [6^,^al^,^6_fc|afc| + b\, ^a}_^,^b-k\aki\ 
+ 92\\{k + k') [4,|6l^,^6_fc|afe| + a\,^b'^_f^,^b_kiaki] 
+ 92±{k + k') [4,^61^, j^6_fe|afc| -h al,^bif,,^b^k^aki] } 

= Y + + al,^bl^,^b-kiaki] 

k,k' 

- 9i±{k k') [4,|6t_fc,|6-fcj,afc| + 4'T^-fc'i^-'=T«fei] 
+ 5-211 (A: + k') [4,^6l^,^6_fc|afct + al,^blf^,^b-kiaki] 
+ 92±{k + k') [4,^6l^,^6_fc|afc| + } 
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k,k' 



t At 



+ g2±{k + A;') [al,^b^_f^,^b-kiak^ + al,^b^_f^,^ b-k^aki] }■ 



(A-4) 



Here we put k' —k' in the second equality. The pairing interaction can be separated into 
spin-symmetric and spin-antisymmetric parts as 

J^{-9i\\{k + k') + g2\\ik + k')} 



H 



mt 



k±' 



+ {-9idk + k')+g2dk + k')} 



^K'T^-fc'i + H'fi-k'^)^{^-kiak} + b-k^aki) 
+ {gi±ik + k') + g2±{k + k')} 



V2' 



The pairing interactions both for triplet and singlet are thus given by 

V^\{k + k') = -g,^^{k + k') + g2ii{k + k'), 
Vl{k + k') = -g,^{k + k') + g2±ik + k'), 
V%k + k') = gMk + k') + g2±{k + k'). 



(A-5) 

(A-6) 
(A-7) 
(A-8) 



Appendix B: Odd-energy superconductivity in Id systems 

In this appendix, wc present a semi-analytic argument to find the most stable type of the 
gap symmetry. Despite its simplicity, it gives the gap symmetry consistent with the numerical 
results in qld system. In order to make our argument as simple and transparent as possible, 
we here discuss the superconductivity in Id, although any kind of phase transition will not 
be realized in purely Id systems due to thermal (and quantum) fluctuations. We assume that 
the pairing interaction have a maximum at q = 2k},^ (see Fig. B-1). Considering that the first 
correction to the vertex have the g-dependence as log \ q — 2kp\ and the renormalized vertex 
is proportional to \q — 2/ci?|°, this assumption is quite natural in Id systems even if we take 
into account the suppressions due to the temperature or the interchain hopping. 

Since the pairing interaction is an even function of momentum with period 27r, we can 
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• »^ ' r^^:^^ jTi • 1 ► 

COS q 

Fig. B-1. Illustration of the pairing interaction V{q) whose maximum is at g = 2kp (solid line). 
The dashed and dash-dotted line indicates cosg and cos2g, respectively. Here we set 2kp — 7r/2, 
assuming the quarter- filling case. 



expand as 



V{q) = ^ + '^i^o^Q + V2COs2q -\ . (B-1) 



Here the coefficients Vn are given by 



71" JO 

They can be written in the separable form 



Vn = — I dqV{q)cosnq. (B-2) 



y*^™^(A:- A:') = ^ +7;icosA;cosA;' + 7;2COs2A;cos2A:' + --- , (B-3) 

Aodd 



V'^'^ik-k') = t;isinfesinA;'-Fu2sin2/csin2A;'-h ••• , (B-4) 

each term corresponding to the pair scattering {k,—k) — > {k',—k'). At quarter filling, the 
coefficients may be estimated as 

vo > 0, ~ 0, < (B-5) 

for the interaction which is peaked at q = 2kF- (The auxiliary cosine curves in Fig. B-1 will 
be helpful for a rough estimate of Vn-'^'^) 

B.l Spin- fluctuation model 

If we adopt the spin-fluctuation model, the pairing interactions in the singlet and triplet 
channels are written as 

~ -I- 3^1 cos A; cos /c' -I- 3^2 cos 2A; cos 2A;' -I- • • • , (B-6) 

~ — f 1 sin /c sin /c' — 772 sin 2/c sin 2A;' -|- • • • . (B-7) 

Prom (B-5), the most attractive component is n = 2 in the spin singlet channel, having the gap 
function cos 2k. If we consider the higher-order components, the pairing interaction peaked at 
q = 2kF will favor singlet pairing with gap function whose general form may be written as 

A"(A;) = ^A„cos2(2n-M)A;, (B-8) 
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which vanishes just at k = kp{= '/r/4), i.e., odd-energy gap function A(A;) oc ^[k) at quarter 
fining. This is because 1)3 0, 1^4 > 0, ~ 0, fg < 0, and so on. In case of half fihing, the gap 
function has the form A(A;) = A„ cos(2n-|- which also vanishes just at A; = kp{= 7r/2). 

The triplet components can also be attractive although subdominant. The higher-order 
component f4 > 0, however, brings moderately large attraction, though V4 < \v2\- Therefore, 
the odd-energy gap function 

A*(A;) ~ ^ A„sin4nA;, (B-9) 

n 

which vanishes just at k = kF{= 7i"/4), is favored also in the triplet channel. We may note 
that as the pairing interaction becomes more singular, becomes comparable to ^2, and 
the difference in between singlet and triplet becomes smaller. For example, the interaction 
V{q) ~ 5{q — 2kF) leads to \v2\ = ^4. 

We have seen that the pairing interaction whose maximum is at g = 2kF generally leads 
to the odd-energy gap function both for singlet and triplet pairing, and the corresponding 
coupling constant for the singlet is (three times or more) larger than that of triplet. If we 
consider the charge fluctuation, this difference is reduced. This general conclusion is completely 
consistent with the numerical results. -"^^ 

B.2 Id Fermi gas model 

Next we consider the Id Fermi gas model with repulsive interaction. From eqs. (A-6)-(A-8), 
the pairing interaction is given by 

V\k + k') = gi{k + k')+g2{k + k'), (B-10) 

V\k + k') = -gi{k + k') + g2{k + k'). (B-11) 

We assume that both gi{q) and (72(9) have a maximum at q = ikp, and expand similarly to 
eq. (B-1) to obtain the separable form 

gi{k + k') = ^^+g^^Koskcosk' + gf\os2kcos2k' + ■■■ (B-12) 
-g^^^ sin k sin k' - gf^ sin 2k sin 2k' , (B-13) 

(0) 

g2{k + k') = ^+g^2^coskcosk' + g^2^cos2kcos2k' + --- (B-14) 
—g^^ sin k sin k' — g^^ sin 2k sin 2k' — • • • , (^^'IS) 



and 



(0) , (0) 

V^{k + k') = +{g^^^ +g'i^) cos k cask' 

+{gf^ + gf^) COS 2k cos 2k' + ■ ■ ■ , (B-16) 
V\k + k') = {g^i^ - g^2^ ) sin k sin k' + (^f ^ - g^^^ ) sin 2k sin 2k' + ■ ■ ■ . 
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(B-17) 

In the same manner as the spin-fluctuation model, we can estimate each component as gf'l > 0, 
^1^2 ^ 9^2 < 0- Therefore, for singlet pairing, the odd-energy gap function is realized as in 
the spin-fluctuation model. For triplet pairing, on the other hand, two types of gap functions 
are possible. First, in the case that the backward scattering is dominant, i.e., \g^^\ > \g^^\, 

(2) (2) 

the attraction originates from n = 2(2m -|- 1) components 5i — ^2 < Oj so that the gap 
function should have the form 

A*(A;) = ^A„sin2(2n + l)A;, (B-18) 

n 

which has a peak at k = kF{= 7r/4), namely, p-wave triplet. Second, in the opposite case, 
i.e., I^i^^l < 152"^!' large attraction comes from n = 4m, and we have the odd-energy gap 
function 

^\k) = ^ AnSin4nA;, (B-19) 

n 

which vanishes at k = kF{= 7r/4). Above results are also consistent with the results of the 
present RG method. The important knowledge here is that the dominant backward scattering 
promote the p-wavc triplet superconductivity, which competes with and can prevail over the 
d-wave singlet one in Id systems. 
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